Optimal parameter estimation of depolarizing channel 
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We investigate strategies for estimating a depolarizing channel for a finite dimensional system. 
Our analysis addresses the double optimization problem of selecting the best input probe state and 
the measurement strategy that minimizes the Bayes cost of a quadratic function. In the qubit case, 
we derive the Bayes optimal strategy for any finite number of input probe particles when bipartite 
entanglement can be formed in the probe particles. 
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I. INTRODUCTION 

In order to design a reliable communication system one 
requires a priori knowledge of the property of a channel. 
Precise knowledge of the channel allows us to devise ap- 
propriate coding, modulation, and filtering schemes. In 
general, the channel property is not stationary, so one 
should first acquire and then track the optimal operating 
point of each device by monitoring the condition of the 
channel. It is important, therefore to know how to esti- 
mate the channel property in an efficient way, that is, as 
precisely as possible with minimum resources. 

A reasonable assumption is that we know that the 
channel belongs to a certain parameterized family, and 
only the values of the parameters are not known. To 
know them one may input a probe system in an appro- 
priate state into the channel and make a measurement 
on the output state. Only when an infinite amount of 
input resource is available, one can determine the chan- 
nel parameters with perfect accuracy. In the quantum 
domain, however, the resource is often restricted for var- 
ious reasons. For example, when one is to monitor a 
fast quantum dymanics at cryogenic temperatures, the 
input probe power should be kept as low as possible so 
as to prevent the system from heating up while obtaining 
meaningful data in a short time. This restricts the avail- 
able amount of probe particles. Furthermore, preparing 
the probe in an appropriate quantum state is usually an 
elaborate process. Thus to find the efficient estimation 
strategy relying only on a restricted amount of input re- 
source is of practical importance. 

In estimating a quantum channel parameter, given a 
finite amount of input resource, both the input probe 
state and the measurement of the output state need to 
be optimized. This double maximization problem has 
been studied in the context of estimation of SU((i) uni- 
tary operation Estimating a noisy quantum chan- 
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nel has been discussed in the literature |[ |j. In 
ref. the locally unbiased estimator and the Cramer- 
Rao bound are extensively discussed for the depolarizing 
channel for a qubit system. The locally optimal strategy, 
which achieves the Cramer-Rao bound at a local point 
of the parameter space was derived when two qubits at 
most are used. This result would be useful in the limit 
of large ensemble of the input probe. In such a limit, 
of course, one can establish the channel parameter with 
a very high degree of accuracy. To improve the rate at 
which the estimation accuracy grows with the number of 
probe particles, one may first apply some preliminary es- 
timation using a part of probe particles to establish the 
most likely value of the parameter, and then use the lo- 
cally optimal strategy around this value to get the final 
estimate ^, |[ 0. Refs. ^ focus on several noisy 
qubit channels. They study some reasonable, although 
not optimal, strategies based on maximum likelyhood es- 
timator, and derive the asymptotic behavior of the cost 
as a function of the number of input probe qubits. 

In contrast, we are concerned here with the Bayes op- 
timal strategy which minimizes the average cost. The 
scenario we have in mind is that one has no particular 
knowledge about the a priori parameter distribution, and 
the available number of probe particles is strictly limited. 
We then take into account the possibility of rather large 
errors. We seek the strategy that works equally well for 
all possible values of the parameter on average, that is, 
the strategy which is more universal for various possible 
situations. 

It seems difficult for us to study this problem for the 
most general probe state. In this paper we deal with the 
depolarizing channel by assuming that we dispose of M 
pairs of probe particles and only bipartite entanglement 
can be formed in each pair. This might be a practically 
sensible assumption from the view point of optical imple- 
mentation given current technology. Our problem is to 
find the best estimation strategy to minimize the average 
cost. We consider the quadratic of a cost function. 
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II. QUBIT CASE 

Let /5 be a density operator in the 2 dimensional Hilbert 
space 7^2- The depolarizing channel Cg maps a density 
operator p to a density operator which is a mixture of p 
and the maximally mixed state, 



Cep = Op 



1 - 



(1) 



The parameter 9 represents the degree of randomization 
of polarization. For the map Cg to be completely posi- 
tive, the parameter 9 must lie in the interval —\<9<1. 

Let us start with two qubit systems as the input probe. 
For simplicity we only consider a pure state family of 
the probe ^ — \^) (^'|. This may be represented in the 
Schmidt decomposition 



I*) = |0) «) |eo> + ^/l-x |1) ® |ei) 



(2) 



where {|0) , |1)} and {|eo),|ei)} are orthonormal basis 
sets for the first and second probe particle, respectively. 
What is the best way to use this state? There are two 
possibilities to consider; 

(a) Input one qubit of the pair into the channel keep- 
ing the other untouched leading to the output state 



■^i{9) = {Ce®i)\^) (*| 



(3) 



(b) Input both qubits into the channel and have the 
output state 



i!2{9) = {Cg®Cg)\^) {^\ 



(4) 



A measurement is described by a probability operator 
measure (POM) t\{9) [|, I), also referred to as a positive 
operator valued measure (POVM) The average cost 
for the quadratic cost function is given by 

C,{x)= I d9 I d9{0~9fz{9)i:r\fi{9)^i{9)\, (5) 



where z[0) is the a priori probability distribution of 0, 
and /_ 1 dOIi{9) = /. It is assumed that we have no a 

priori knowledge about 9, that is, z{9) ~ |. Given the 
channel £9, we are to find the optimal probe \^) and the 
POM ti.{9) minimizing the average cost C{x). 
It is convenient to introduce the risk operator 

W{9) = - [\d9'{9-9'fi>,{9'), (6) 
= _ 20^(1) +6)2^(0)^ (7) 

where W^^^ = IJ^ 1 d99^'^i{9). The average cost is then 



C{x) = TrF, r = / d9t{{9)W{9) 



(8) 



For a fixed probe state |*), the optimal POM 11(6*) is 
derived from the necessary and sufficient conditions to 
minimize the average cost |l^ : 



(i) f = ft,and 



W{9) - f tl{9) = for aU 



(ii) W{9) - f > for aU 9. 

The optimal solution for a single parameter estimation 
with a quadratic cost is well known |l^ . The optimal 
POM is constructed by finding the eigenstate \9) of the 
minimizing operator O which is defined by 



(9) 



that is, n{9) = \9) {9\ so that e\9) =9 \9). We then have 
f = _ 91^(0)0 fi-ojj^ ^jjjgjj the conditions (i) and 

(ii) are easily verified. 

For a discrete system, one can find the optimal POM 
with finite elements. Let the spectral decomposition of 
^(0) fQj. Q^j. two-qubit system be 



Then the minimizing operator is 
2 , , 



e 



4 



(10) 



(11) 



Let the spectral decomposition of 8 be 

e 

The optimal POM is then given by 

4 

n{9)=Y,S{9-9,)\9,} {9, 

i=l 



4 
i=l 



(12) 



(13) 



This implies that the measurement has 4 outputs at most 
and we then estimate the channel parameter as one of 4 
9i's. Before going on to derive the optimal strategies, 
let us define some notations. As seen below the output 
states ^i{9)^s can be written as a direct sum 



(14) 



where V'i(^) is in the subspace spanned by \p,i) = 
|0) 1/0) and 1^2) = |1) 1/1), and (j>i{9) in the subspace 
spanned by = |0) (g) |/i) and I1/2) = |1) «> |/o). In 
the following all 2 x 2 matrices represent density operators 

in H^, with \pi) = ^ q ^ and \p2} = 1 

Case (a): 

The output state 'i'i{0) is given by 



M0) 



il + 0)x 29y/x{l-x) 



(15) 



29y/x{l - x) (l + 6l)(l-x) 
^\l^x)\i^i){,yi\+x\i^2){'y2\ ■ (16) 
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The elements of the risk operator are 
2x 



27 V 



iy(2) 



27 
1 

27 



- x) 

^x{l-x) 2(1 -x) 

8x 7^a;(l - a;) 
_ 7^Jx{l-x) 8(1 -x) 

6x 5^a;(l — a;) 



b^Jx{l-x) 6(1 -x) 



©^ij, (17) 
®^i)(18) 
®^i)(19) 



where (fi — [1 — x) \vi) + x \v2) {i^2\- After a lengthy 
but straightforward calculation (see Appendix ^) we 
have 



l+X 2yjx{l - X) 

2^x{l - x) 2-x 



1 - 

9^^ 



(20) 



To diagonalize we introduce r = ^Jl + \2x(l^^x) and 

(21) 



C0S7 



r - 1 + 2x 



sin7 



r + 1 - 2a; 



2r ' V 2r 

The eigenstates and eigenvalues are then 

\ei) = C0S7 l^i) + sin7 1/^2) , 6'i = (3 + r)/9, 

16*2) = -sin7 l^i) + C0S7 1^2) , 6*2 = (3 - r)/9, 

1^3) =|^^i), 03 = 1/9, 

\0i) =|J^2), 04 = 1/9. 



The average cost finally reads 

Ci{x) = Tr(W^(2) _ qw'^o)q^ ^ 8 

81 



(22) 



l + (-^)^ 



. (23) 

This is minimized by the maximally entangled state input 

|vI/) = -^(|0)®|/o) + |l)®|/i)), (24) 

for which = | and O2 ~ 0^ = 0^ = ^. Therefore the 
optimal measurement is actually constructed by the two 
projectors 



ni = |'j')(*|, n2 = /-|*)(*| 



(25) 



with the associated guesses = | and 6*2 = ^, respec- 
tively. The minimum average cost is C'lmin = 
Case (b): 

The output state ^2(0) = "02 (0) © 4>2(0) is given by 



- {\ -x)e + e'^ e'^^x{i-x) 
e-^^x{i-x) \ + {i-x)e + e-' 



1 



MO) = 
MO) = 

The elements of the risk operator are 



(26) 
(27) 



1 

27 



4 + 9a: 7^x{l-x) 
7^a;(l - a;) 13 - 9x 



5 . 

27^*' 



1 

27 
1 

27 
1 



7x b\/x{\ — x) 
b^x{l-x) 7(1 -a;) 



4 + 75x &l.Jx{l~x) 

405 [6lVa;(l"a;) 79 75a: 
11 . 
® 405^^- 

The minimizing operator is (see Appendix ??) 



e = 



17[13 + 8a:(l - a;) 



a c 
c b 



(28) 



(29) 



(30) 



(31) 



where 



a = 7a;(35-20a; + 2x2) 
b = 7a;(17+16a; + 2a;2) 



c = 9[9-2a;(l-a;)]^a;(l-a;). 



(32) 



To diagonalize it we use r = y^a — 5p + 4c? and 



COS7 



r + a — b 



2r 



- , sin7 



r — a + 6 



2r 



(33) 



We then have the similar eigenstates to Eq. (E2h and the 
^^^^TitroiTi^ri n — a n — a n — a — t with 



eigenvalues 9i = 9+, 62 = 0_, and 



y± = 



119[l + 2x(l-a;)]±; 



34[13 + 8a:(l - x)] 
The average cost is then 



(34) 



8[391 + 606a:(l - x) - 10a;^(l - a:)^] 
^""^ - 2295[13 + 8a;(l-a;)] ' ^^^^ 

This is an upward convex function, symmetric with re- 
spect to a: = ^. The minimum is attained at x = 0, 1, 
that is, by separable input states. This reads (72111111 = 

184 
1755 ■ 

The average costs for cases (a) and (b) are shown in 
Fig. |l|: 6*1(3:) (sohd hue) and C2ix) (dashed hne). We 
see that C'lmin < 62111111 so that the optimal estimation 
strategy, using two probe qubits, is to prepare them as a 
maximally entangled pair and to input one qubit of the 
pair into the channel keeping the other untouched. The 
estimation is then obtained by applying the two element 
POM, Eq. (^, as described in Case (a). This strategy 
is represented schematically in Fig. |[ 

When M maximally entangled pairs j^i)®*^ are avail- 
able, it is best to use them so as to have the output 
[{Ce (g) /) 1^') (\1>|]®*^. The optimal measurement for this 
can be derived straightforwardly. This is discussed in the 
next section as a part of an arbitrary finite dimensional 
case. 



4 



-4— ' 


0.13 


o 




o 


0.12 








0.11 








0.10 















\ / 
\ / 

/ \ 






>/ 


/ \ 






—/^ — \ — 











0.0 



0.25 0.5 

X 



0.75 1.0 



|M/)< 



FIG. 1: The average costs as a function of x 




□ i = |vl/)(vl/| 

n2 = ;-|vi/)(M/| 



9 = 1 
" = 1 



FIG. 2: The optimal estimation strategy using two probe 
qubits. I^*) is the maximally entangled state. The output 
state is projected onto {111,112}. We guess the channel pa- 
rameter as 6 — ^ for the outcome IIi and ^ = | otherwise. 



III. d-DIMENSIONAL CASE 

The action of the depolarizing channel on a d dimen- 
sional system is described by 



Cgp = 9p 



1 



-I. 



(36) 



Complete positivity then implies — ^p^j < < 1. For 
d > 3, we have not succeeded in finding the optimal probe 
state, even when we restrict ourselves to a pure state. In 
this section we focus on the most plausible input state, 
that is, the maximally entangled state, and consider the 
estimation using M entangled pairs. Only for d — 2, is 
the optimality ensured. 

It might be interesting to compare the three cases spec- 
ified by the three different outputs; 

(a) M product states of the pair 

^i{e) = {Ce(E>i)\^){-^\, 



where I^E") is the maximally entangled state, 
(b) M product states of the pair 



(37) 



1 _ fl2 



d2 



-I^I, 



(38) 



(c) 2M product states of 

m = Cg\0) 



(39) 



(The input state in case (c) can be any pure state in the 
d dimensional space.) Let us first consider the case (a). 
We denote Eq. iM) as 



where 



and 



*(0) = /o(0)do + /i(e)ai, 



ao = |*)(^'|, ai = /-|'J')(*| 



fo{e) = 



1 



hid) 



1 



d2 ' ' d2 

The output state can then be represented as 

M 

^(^)«5M ^ ^ /o(0)^^-™/i(0)"i„ 



m— 



where 



Ara= ^ Oil (g) • • • (g) ai^^ , 

(n + ...iA/=m) 



(40) 
(41) 

(42) 

(43) 
(44) 



is the projector onto the symmetric subspace. The risk 
operator is 



M 



W{6) ^ t^™' - + (45) 



m=0 



where w};:^ ee J\ dOO'' fo{e)^^^''' fi{9)"^ ■ The optimal 
POM is 



M 



m—O 



where 



(46) 



(47) 



We then note that 

W{0) - f - ^ (0 - 0„O'w^n^^m > 0, (48) 



M 



from which it can easily be seen that the conditions (i) 
and (ii) hold. The minimum average cost is 



M 



m=0 



,(2) 



(0) 



(d2_l)". (49) 
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The other cases can be dealt with in a similar manner. 
In the case (b), we just put 



0.15- 



fo{9) = 



1 



d2 



^2 1 

h{0) = - 



cP 



(50) 



The minimum average cost C2{M) is then given by the 

same expression as Eq. (|9|) with w^^'s defined by fo{e) 
and h{e) of Eq. @. 
In the case (c), we use 



ao 



ai = /-|0) (0|, 



and 



1 - 1 



1 - 



(51) 



(52) 



The minimum average cost is 



C'sEp(M) 



2M 

E 

m— 



,(2) 



2M 
m 



[d-iy 



(53) 

The three costs Ci(Af), C2(M) and Csep(M) are 
plotted in Fig. | (d = 2), Fig. | (d = 3), and Fig. 
^ (d = 10). In the figures another average cost Cmi.{M) 
is also plotted. This cost is by the strategy belonging to 
the case (c), but unlike the one attaining C'sep(M), the 
estimator is made by the maximum likelyhood principle 
for which 



md 



(54) 



instead of Eq. (^ 

C'ml(M) = 



2M{d - 1) ' 
and leads to the analytic expression 
1 d'^id + i) 



2M 6(d2 - 1)3 



(55) 



It is this strategy that was used in ref. Q] for the case of 
d = 2. 

For d > 3, the minimum average cost is always attained 
by a separable probe state. Only in the two dimensional 
case, is it the bipartite entangled probe that attains the 
minimum average cost. It is worth mentioning the de- 
polarizing channel with the narrower parameter region 
< < 1, which is a more commonly used model with 
an well defined interpretation of randomized probability 
of 9. We found that the best probe in this model is al- 
ways a separable state. In this sense a separable state 
is generally an adequate probe state for the depolariz- 
ing channel estimation as far as the comparison with a 
bipartite entangled probe state is concerned. 
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FIG. 3: The average costs as a function of the number of 
pairs. 
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FIG. 4: The average costs as a function of the number of 
pairs. 
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IV. CONCLUDING REMARK 

When we have several identical samples at our dis- 
posal, it might be desirable to apply the best collective 



FIG. 5: The average costs as a function of the number of 
pairs. 
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measurement on the whole system. This means prepar- 
ing a single multi-qubit state followed by an optimized 
measurement. We might also consider performing a pre- 
liminary measurement on a part of the system and then 
fccdback this back to deal with the remaining part. But 
in the case of the previous section, the collective measure- 
ment on M identical output pairs or 2M identical output 
particles is not necessary. The action of the depolarizing 
channel on a maximally entangled state always results in 
a statistical mixture between the input state and its or- 
thogonal complement (Eq. (|40|)). Estimating the channel 
parameter is nothing but determining this mixing ratio, 
which is a classical distribution. Therefore the optimal 
measurement is realized by a separable type constructed 
by the binary orthogonal projectors {ao,ai} according 
to Eq. (^). In the case where the output state includes 
the channel parameter as a quantum distribution, that 
is, the parameter appears in the off diagonal components 
in the density matrix, the optimal measurement would 
be a collective measurement. When the channel includes 
a unitary opreration, we will have to face this problem. 
Channel estimation for such a case is a future problem. 

It is a remaining problem to see how effective the mul- 
tipartite entangled probe is. However, in the estimation 
of decoherence channel under the power constraint sce- 
nario, that is, under a given and fixed number of probe 
particles, it seems more common that entanglement is not 
necessary. In fact, in the cases of the amplitude damp- 
ing channel and dephasing channel, there is no merit to 
use entangled probe. In the amplitude damping channel, 
for example, the best probe is to input the most highly 
excited state. An entangled probe is rather wasteful be- 
cause this includes the state components other than the 
excited state and these components are less sensitive to 
the damping. 

Finally it might be interesting to study the multi pa- 
rameter case, such as the Pauli channel estimation. We 
may then ask how to optimaize (in Bayesian sense) the 
simultaneous measurement on the noncommuting observ- 
ables as well as searching for appropriate probe states. 



APPENDIX A: DERIVATION OF EQ. (20) 



For obtaining the minimizing operator Q in Eq. (|2C|), 
we first diagonalize iy(°) by Uq — uq® where 



COS70 -sm7o 
sin7o COS70 



ith r-Q = — 3.t(1 — x) and 



COS70 



2r, 



sm7o 







ro + 1 - 2a; 



2rn 



The spectral decomposition 

4 



(Al) 



(A2) 



(A3) 



is given by 



ki) ^woImi), wi = (l + ro)/3, 

1^2) = Uq |/i2) , W2 = (1 - 7'o)/3, 

ks) = , (1 - .x)/3, 

1^4) = \V2) , = x/'d. 



(A4) 



We then calculate 

4 



(A5) 



where 



This gives 



where 



|1> 

|2) =|M2), 

|3) 

|4) =^2). 



y 



(A6) 



(A7) 



4ro(l+ro)-|-3a:(l-a;) S,(l-2x)^J x{l-x) 

ro(l+ro) ro 

?,(l-2x)^x{l-x) ira(l~ro)-3x(l-x) 

ro ro(l~ro) 



(A8) 
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APPENDIX B: DERIVATION OF EQ. (|^) 

The unitary operator for diagonalizing W^^"^ in Eq. 
^) is Uo — iiQiS I4, where 



Mo 



COS70 — sm7o 
sin7Q COS70 



(Bl) 



7 



with To = y/81 - 128x(l - a;) and 



We then have 



/ro-9(l-2a;) . /ro + 9(l-2a;) 

COS70 = \ , sin7o = W - 



2ro 



The spectral decomposition 



2ro 







4 








1=1 




is given by 








ki> 


= "0 Ia^i) 


LJl = 


(17 + ro)/54, 


1^2) 


= "0 \^J■2) 


, = 


(17-ro)/54, 


ks) 


= Wi) , 


CJ3 = 


5/27, 


1^4) 


= W2) , 


UJ4 = 


5/27. 



(B2) 
(B3) 



6 = 9^,® 



where 



7[ro+9-16a:(l-x)] 8(l-2x)^a;(l-a:) 

ii(l-2x)^x{l-x) 7[ro-9+16z(l-z)] 
17ro ro(17-ro) 



(B5) 



(B6) 



(B4) Substituting this to 6 = U^QUl, we have Eq. (|l]). 
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